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Abstract 

The dipole flexoelectric (FE) polarization in liquid crystals is de- 
rived in the thermodynamic limit at small distortions and small den- 
sity. General microscopic expressions for the FE coefficients are ob- 
tained in the case of the uniaxial and biaxial nematic phases composed 
of C2v molecules. The expressions involve the one-particle distribution 
function and the potential energy of two-body short-range interac- 
tions. In the case of the biaxial nematic phase, six basic deformations 
produce FE polarization but there are only five independent FE coef- 
ficients. 

PACS number(s): 61.30.Cz, 77.84.Nh 

Biaxial nematic liquid crystals are characterized by anisotropic positional 
short-range order and orientational long-range order [1]. The anisotropic 
molecules tend to be parallel to selected axes, labeled by the unit orthogo- 
nal vectors L, M, and A^. In the uniaxial nematic phase only the N axis is 
defined. Stable biaxial phases were observed in 1980 in lyotropic systems by 
Yu and Saupe [2] but their existence in thermotropic systems was not certain 
for many years [3]. Recently Madsen et al. |1] synthesized banana-shaped 
mesogens and evidence for biaxiality was achieved using NMP spectroscopy. 
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Acharya et al. [5] revealed biaxiality of bent-core mesogens by means of low- 
angle x-ray diffraction. Merkel et al. [B] carried out infrared absorbance mea- 
surements on two liquid crystalline organo-siloxane tetrapodes and showed 
the existence of a biaxial nematic phase. 

Theoretical predictions of biaxial nematic phases started in 1970 with the 
paper by Preiser [7]. Then they were studied using a number of theoretical 
methods, e.g. mean- field theory [H], [S], counting methods [lU], [TT], Landau- 
de Gennes theory [12], [13], bifurcation analysis [H], and density-functional 
theory [15]. All the theories mentioned above predict that the system will 
exhibit four phases as the molecular biaxiality varies: a positive and a neg- 
ative uniaxial phases, respectively, formed by prolate and oblate molecules, 
a biaxial, and an isotropic phase. The nwmatic-isotropic phase transition 
is expected to be first order and to weaken as the biaxiality increases until 
it becomes continuous at the point (Landau bicritical point) of maximum 
molecular biaxiality. At this point the system should go directly from a bi- 
axial to an isotropic phase. The uniaxial-biaxial transition is expected to be 
second order. 

A microscopic mean field theory [TB] , [IT] predicts the possibility of lines 
of second-order and first-order uniaxial-biaxial transitions joined at a tricrit- 
ical point. The experimental results by Merkel et al. [6] were successfully 
interpreted in terms of this theory. A weakly first-order uniaxial-biaxial tran- 
sition was also revealed by dynamic light scattering from orientational order 
fluctuations in a liquid crystalline tetrapode [TB]. Recently, two-particle clus- 
ter theory was applied to study the biaxial molecules in the Sonnet model ^16] 
and qualitative agreement with the mean-field theory was obtained [19]. A 
Monte Carlo study of biaxial nematic phases composed of V-shaped molecules 
was done by Bates and Luckhurst [2D]. 

In most practical circumstances the liquid crystal phase alignment is de- 
formed. The deformations usually are described by a continuum theory where 
the free energy is expanded to the second order around the free energy of 
the undeformed state in ascending powers of a parameter that measures the 
deformation. The free energy due to the distortion of the axes is expressed 
in terms of the vector spatial derivatives and the elastic constants. 

In a deformed uniaxial nematic liquid crystal, there should appear in 
many cases a spontaneous dielectric polarization described by Meyer [2T] : 

P = eiiV(V ■ N) + csiN ■ V)N = eiiV(V ■ N) - csN x (V x iV), (1) 

where Ci and 63 are the splay and the bend fiexoelectric (FE) coefficients. 
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respectively. The appearance of spontaneous polarization in liquid crystals 
as a result of orientational deformations is referred to as the flexoelectric 
effect. 

A microscopic mechanism of the FE effect was proposed by Meyer |2T] . 
who pointed out that under the condition of orientation deformation of a 
liquid crystal, banana-shaped or conical molecules should be so oriented 
that their constant dipoles are ordered and macroscopic polarization sets 
in. On the other hand, Prost and Marcerou [22] showed that polarization 
in a deformed hquid crystal is also produced as a result of a gradient in the 
average density of the molecule quadrupole moments. Such a polarization 
does not need asymmetry in molecular shape of the banana or cone kind. 
Later, the FE coefficients for the uniaxial nematic phase were calculated by 
means of a statistical- mechanical theory [23], mean-field theories [2l], [25], 
the density-functional formalism [26], [27], an integral equation approach 
[28], and computer simulations [29], [30]. The expressions connecting the 
molecular asymmetry, the elasticity of the liquid crystal, and the FE coeffi- 
cients were obtained by Derzhanski and co-workers [21]. It was also shown 
[52] that the FE coefficients can have nontrivial dependence on the details of 
the molecule's chemical structure (an isomerization, a charge distribution). 

The flexoelectric polarization can influence electrooptical properties, de- 
fect formation, and structural instability. Thus, different techniques have 
been suggested to observe possible mechanisms producing the FE effect [33] • 
The FE coefficients can be obtained experimentally from measuring the po- 
larizations or the surface charges induced by an imposed distortion or using 
the inverse effect, because when an electric field is applied on a nematic, the 
alignment may become distorted and this will imply a polarization [31] . The 
flexoelectric polarization of 5CB was measured by means of a pyroelectric- 
effect-based technique [3S], [3S] and recently a technique inspired by the 
flexoelectric-optic effect was demonstrated [37] . 

Below, we would like to investigate the FE effect in biaxial nematic 
phases. In the case of the uniaxial nematic phase there are three independent 
phase deformations: splay, twist, and bend. The symmetry considerations 
of Rudquist and Lagerwall [3H] prove that in uniaxial nematic liquid crystals 
the FE effect can be induced by splay or bend. The polarization connected 
with bend has to be perpendicular to the director A^, whereas in the splay 
deformation, a polarization along the director is admitted. The twist is not 
connected with a local polarization of the medium because there is always a 
two-fold symmetry axis perpendicular to the helix axis. 
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It was shown [22] that in the case of the biaxial nematic phase there are 12 
independent phase deformations and six of them are connected with splays 
and bends of the vectors L, M, and N. Thus, we can expect that there are 
six FE coefficients for the biaxial nematic phase. 

Let us assume that the FE polarization, which is a vector quantity, de- 
pends on the spatial derivatives of the vectors L, M, and N 

Pa = SijRiadpRjfS + bijRipdpRja, (2) 

where we denoted 

La = Rla^ Ma = R2a^ = R^a- (3) 

The matrix elements Ria (i = 1, 2, 3 and a = x,y, z) satisfy the conditions 

RjaRja = ^ij) RigRiff = ^a0- (4) 

a i 

The above relations express the orthogonality and the completeness of the 
local frame. It is also possible to derive the identity (to be used later) 

^{RiadpRip + RipdpRia) = 0. (5) 

i 

The number of independent FE coefficients will be determined with the 
help of symmetry requirements. The local frame can be transformed into the 
new one 

R'ia = TijRja, (6) 

where Tij = 1,2,3) are the elements of the orthogonal transforma- 
tion. The polarization may be expressed in terms of new variables with 
new (primed) FE coefficients. As the FE coefficients do not change under 
symmetry operations, we can identify the independent FE coefficients. As a 
result we get six FE coefficients sa and bu. But from the identity ([5]) we con- 
clude that the coefficients su and bu are not unique. The same polarization 
P can be obtained by means of the transformed coefficients 

4i = Su + c, = bu + c, (7) 

where c is any constant. Thus, only five coefficients are independent. Note 
that the differences Su — Sjj, bu — bjj, or Su — bjj do not depend on the constant 
c. 
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Now we consider six small deformations of the directors (L, M, N) of the 
form 02] 

L{r} = [l.q'jy - q5X,qgZ - q^x], 
M{r) = [q^x - q7y,l,qsz - qay], 

N(r) = [q^x - q9Z,qQy - qsz,l], (8) 

where the parameters (i = 4, . . . , 9) describe the deformations {l/qi is a 
certain length much greater than the size of the sample). The corresponding 
FE polarization has the form 

Px = qrisii - 622) + qgisu - hs) = Qio-i + ^gOg, 

Pz = 54(533 - bu) + g6(s33 - ^22) = 9404 + geoe, (9) 

where we introduced six physical FE coefficients aj which satisfy the identity 

a4 — — + a-j + ag — ag = 0. (10) 

Eqs. are phenomenological expressions and we should provide alternative 
microscopic expressions in order to obtain microscopic expressions for the FE 
coefficients. 

Let us consider a system of N rigid molecules with the C2v symme- 
try. Such molecules can form uniaxial {Dooh) and biaxial {D2h) nematic 
phases. The free energy for the system can be derived in the thermodynamic 
limit (A^ oo, V — >■ oo, N/V = const) from the Born-Bogoliubov-Green- 
Kirkwood-Yvon hierarchy [lOj or as a cluster expansion [H]. The total free 
energy F consists of the entropy term and the interaction term 

P = Pent + Pint, (H) 

where 

PPent = J dfdRGif,R){\n[Gif,R)A]-l}, (12) 

PPrnt = -^J df,dR,df2dR2Gin,R,)Gif2,R2)fl2. (13) 

We define dR = dcpdO sin 9dip , /12 = exp(— /?$i2) — 1 (the Mayer function), 
P = l/{kBT), u = r2 — ri = uA; and A is related to the ideal gas properties. 
The normalization of the one-particle distribution function G is 

J drdRGif, R) = N. (14) 
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The equilibrium distribution G minimizing the free energy fllip satisfies the 
equation 

ln[G'(fi, Ri)K] - J df2dR2G{f2, i?2)/i2 = const. (15) 

In the homogeneous biaxial nematic phase composed of or molecules 
the distribution function has the form [39] 

Go{R) = Go{l- L, I- N,n-L,n-N), (16) 

where the unit orthogonal vectors (/, m, n) describe the molecule's orienta- 
tion. For the C2V molecules, the molecule symmetry axis is determined by 
the vector n. Note that this is the long axis of wedge-shaped molecules and 
the short axis of banana-shaped molecules. In order to derive expressions for 
the elastic constants it is enough to assume that, in the deformed phase, the 
phase orientation depends on the position but the magnitude of the alignment 
is constant, 

Go{f,R) =Go[l- Hf^J- N{f),n- L{f),n- N{r}]. (17) 

In order to derive expressions for the FE coefficients we have to take into 
account a small change of the alignment, 

G{f,R)=Go{f,R)[l + g{r,R)], (18) 

where g is expected to be small. The microscopic polarization depends on 
the position inside the phase via the distribution function 

P{f) = J dRG{f,R)jl{R). (19) 

The molecule electric dipole moment can be defined in the molecular frame 
as 

l^a{R) = f^iL + /i2"^a + fisria, a = x,y, z. (20) 
According to Straley [23], the function g can be obtained from the expression 

g{f, R,) = I dudR2h2{u ■ V)Go(f, R2). (21) 

Finally, from Eqs. (fT9|) . (pTI) . and ([8]), we get the components of the micro- 
scopic FE polarization, 

Px = J dudRidR2fl2GoiRl)fJ-xiRl)[iU2z-W2x)q9Uz + U2yq7Uy], 

Py = J dudRidR2f 12GoiRl) fiyiRl)[W2yi-q8Uz) + U2yi-q5Ux)], 

= J dudRidR2fl2Go{Rl)lXz{Rl)[iU2z-W2x){-q4Ux) + W2yq6Uyl;22) 
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where we denoted 

Ua = diGola + dsGoUa, Wa = 92G'o^a + d^GoUa. (23) 

When we compare Eqs. and fl22p . we obtain the equations for the FE 
coefficients: 

04 = J dudRidR2fi2Go{Ri)i2,{Ri){-u^){U2z-W2x), 

= J dudRidR2fi2Go{Ri)ny{Ri){-Ux)U2y, 

ag = J dudRidR2fi2GoiRi)Hz{Ri)uyW2y, 

a7 = J dudRidR2fi2Go{Ri)n:,{Ri)uyU2y, 

as = J dudRidR2h2Go{Ri)iiy{Ri){-Uz)W2y, 

ag = 1 dudR,dR2fi2Go{Ri)f^ARi)uz{U2z - W2.). (24) 

Let us consider the uniaxial nematic phase composed of wedge-shaped molecules. 
The long molecule axes n are almost parallel to the axis, and we get 

= 0, as = a7 = 0. (25) 

We recover the Meyer expression ([T]) with two FE coefficients, 

ei = a4 = ag, 63 = -ag = -ag. (26) 

Finally we consider the uniaxial nematic phase composed of banana-shaped 
molecules. The long molecule axes / are almost parallel to the L axis, and 
we get 

Wa = 0, as = ag = 0. (27) 

We obtain the Meyer expression ([T]), where is replaced with L. The FE 
coefficients are 

ei = ay = ag, 63 = — a4 = —a^. (28) 

Note that Eqs. fl26|) and fl28l) describe C2V molecules, whereas very often 
simpler expressions for Coov molecules are present in the literature. 

In conclusion, we derived the microscopic formulas for the six FE coef- 
ficients in the case of the biaxial nematic phase with the C2V molecules. It 
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appeares that only five FE coefficients are independent. In order to calcu- 
late the values of the FE coefficients one needs the one-particle distribution 
function and the potential energy of molecular interactions. The Meyer ex- 
pressions [21] are recovered in the case of the wedge-shaped and banana- 
shaped molecules in the uniaxial nematic phase. Generally, the splitting of 
the two Meyer FE coefficients and the appearance of new small FE coeffi- 
cients are expected at the uniaxial-biaxial nematic transition. In order to 
describe real substances, the presented results should be generalized beyond 
the low-density limit, where the Mayer function /12 is replaced with a bet- 
ter approximation of the direct correlation function C2. On the other hand, 
other sources of dielectric polarization should be taken into account: the 
quadrupole contribution or the ordoelectric polarization. But even then, a 
qualitative comparison between theory and experiment may be difficult, be- 
cause the experimental data on FE coefficients are still scarce and sometimes 
contradictory [33] . 
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